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Wo quantify the extractablc entanglement of excited states of a Licb-Linigcr gas that are obtained 
from coarse-grained measurements on the ground state in which the boson number in one of two 
complementary contiguous partitions of the gas is determined. Numerically exact results obtained 
from the coordinate Bethe ansatz show that the von Neumann entropy of the resulting bipartite 
pure state increases monotonically with the strength of repulsive interactions and saturates to the 
impenetrable boson limiting value. We also present evidence indicating that the largest amount of 
entanglement can be extracted from the most probable projected state having half the number of 
bosons in a given partition. Our study points to a fundamental difference between the nature of the 
entanglement in free-bosonic and free-fermionic systems, with the entanglement in the former being 
zero after projection, while that in the latter (corresponding to the impenetrable boson limit) being 
non-zero. 
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I. INTRODUCTION 

In the quantum information and communication com- 
munity, studies of the properties of quantum entangle- 
ment and entanglement measures are mainly motivated 
by the idea that entanglement is a vital resource for pro- 
cessing quantum information. Of particular current in- 
terest is the quantification of the entanglement present in 
many-body quantum systems [1] for which experimental 
measurement protocols have been recently proposed [2]. 
In this paper we will study the entanglement that can be 
extracted from bipartite projections of the ground state 
of the integrable repulsive Lieb-Liniger model on a ring 
[3, 4]. The corresponding hamiltonian describes a con- 
tinuous one-dimensional gas of bosons with short-range 
interactions modeled by Dirac delta- functions. Though 
these one-dimensional systems have long been studied in 
atomic waveguides with tight axial confinement [5-7] , the 
ring-shaped traps needed to reproduce periodic boundary 
conditions have only been recently developed [8]. 

The main objective of our work here is to study the ef- 
fects of the strength of interactions on the entanglement 
obtained from projections of the ground state of this con- 
tinuous many-body system — an aim that is facilitated by 
the fact that the Lieb-Liniger hamiltonian is exactly solv- 
able by the Bethe ansatz. In fact, the entanglement en- 
tropy in other continuous integrable systems, for example 
the Calogero-Sntherland model [9, 10] and the anyonic 
Lieb-Liniger model [9, 11, 12], has already been studied 



under the framework of particle partitioning, but as we 
explain in the following paragraphs our extraction proce- 
dure involves a spatial partitioning of subsystems. More 
recently, the von Neumann and Renyi block entropies in 
the ground state of spatially continuous free-fermionic 
models [13] were obtained from conformal field theory 
arguments as the appropriate limit of the corresponding 
lattice model results [14, 15]. 

Although there is as yet no conventional approach to 
quantifying the amount of entanglement in multi-particle 
assemblies of indistinguishable particles [16-18], the type 
of projections we consider here allows us to give a definite 
measure of entanglement in the resulting projected states 
via the von Neumann entropy. In such indistinguishable 
particle systems, entanglement between two spatial par- 
titions of the system may both arise from the correlation 
in occupancies of the partitions, as well as the corre- 
lation between distributions of the particles in the two 
partitions. For example, occupancies of two boxes, by a 
total number of two indistinguishable bosons can be of 
the form ]2, 0) + Jl, 1) -I- JO, 2),^ and there is entanglement 
in this state by virtue of the uncertainty of the number of 
particles in each box. If one can somehow eliminate this 
form of (somewhat trivial) entanglement, then does one 
still retain any entanglement between partitions in a sys- 
tem of indistinguishable particles? We will show below 
that the answer to this question is "no" for free bosons, 
while it is "yes" for free fermions. providing a striking 
qualitative difference in the cause of the entanglement in 
these systems. 
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In this paper, wc arc interested in quantifying the en- 
tanglement in projections of the Lieb-Linger ground state 
that are obtained after coarse-grained measurements re- 
veal the number of particles in one of two contiguous par- 
titions A and B that divide the gas into two half-rings. 
Thus, the projected pure state is spatially partitioned 
(bipartite) and one can quantify the resulting entangle- 
ment by the von Neumann entropy. A similar procedure 
has been used to quantify the entanglement extractable 
from stationary and non-stationary states of impenetra- 
ble boson gases [19] , supersinglet states, and several spin 
chains [20] . In the spirit of these latter works we study 
here the notion of a projectively extractable pure state 
entanglement £pp as defined by 



max {p{k)SA{k)} = maxjffe}, 

k k 



(1) 



where p{k) is the probability of projecting the ground 

state into a state XAB{k) having the fraction k/N of par- 
ticles in region A and SA{k) = ~ tr[pA{k)\og2 pA{k)] is 
the von Neumann entropy of the reduced density matrix 
PA{k) = trs \XABik)){xABik)\. That is, if 11^ is a projec- 
tor onto the state subspace having k particles in partition 
A, then the projected ket we are interested in may be ex- 
pressed as \xAB{k)) = n^lx) in terms of the groimd state 
ket Ix). The entanglement measure above (1) gives the 
maximum weighted entanglement £k = p{k)SAik) over 
all possible projection outcomes and consequently cap- 
tures the probabilistic nature of the preparatory mea- 
surements. 

Our analysis begins in Sect. II with a brief introduc- 
tion to the coordinate Bethe ansatz that forms the basis 
of our exact computations. The probability of obtaining 
each projection outcome is calculated in Sect. Ill, while 
the von Neumann entropy of the projected states and 
their corresponding weighted entanglement are given in 
Sect. IV. On the basis of these results we argue that the 
projectively extractable pure state entanglement £pp is 
equal to the weighted entanglement £n/2 of the balanced 
case in which exactly half of the bosons are present in 
both partitions (the total number TV is given to be even) . 
Our main conclusion is that the extractable entanglement 
£pp increases monotonically with interaction strength 
and saturates to its impenetrable boson value. We con- 
clude with a summary and some remarks in Sect. V. 



II. BETHE ANSATZ 

We briefly review here some of the well-established 
properties of the Bethe ansatz solution for the ground 
state of the Lieb-Liniger model [3, 21]. Written in 
dimensionless form, in which lengths are measured in 
units of the ring circumference L and energy in natural 
units h'^/2mL'^, the Schrodinger equation for N delta- 
interacting bosons is 



^x(x) 



N 

j=l i l<k<j<N 



X(x). (2) 



The dimensionless interaction constant is taken to be 
non-negative c > so that the gas docs not collapse. 
The unique periodic ground state eigensolution of this 
hamiltonian is given by the normalized coordinate Bethe 
ansatz 



{V} 



(3) 



Fp(x) 



n 



X-Pj - Xpk - icsgn{xj - Xk)_ . . 

, 1/2 ' W 



<k<j<N {^!nfc<,[(A,-A,)2+c2]} 



where the momentum vector Xp has A'^ components Xpj 
that form a permutation V of the solutions Xj of the 
Bethe equations 



jV 



n Aj - Xk+ ic 
Xi - Afc - 



fe=i 



tc 



(5) 



which satisfy Aj = 0. The quantity (— 1)^^' is the 
signature of the permutation V and the normalization 

factor Af may be obtained by, for example, the quantum 
inverse scattering method [22]. Explicitly, the absolute 
square \J\f\ is the determinant of the second derivatives 
of the Yang action S evaluated at the solutions of the 
Bethe equations (5): 



det 



dXjdXk 

N 



= det \^Sjk + ^ 



2cS 



'jk 



~i (Afc - XiY + c2 
2c 

~ (A, -Afe)2+c2 



(6) 



III. PROJECTION PROBABILITIES 

In this section we are interested in calculating the prob- 
ability p{k) of successfully projecting the ground state 
X into the pure state XAB{k) having k bosons in par- 
tition A and N — k bosons in partition B. For conve- 
nience let us define the vectors = [xi,. . . ,XkY and 
xs = {xk+iT ■ ■ tXnY . Since the wavefunction x(x) is 
manifestly invariant over any permutation of coordinate 
indices, we can express the success probability as 



p{k) 



I I lx(x) 

J A J B 



' dx.Bds.A- 



(7) 



Also, because of the translational invariance of the 
ground state wavefunction we may choose the partitions 

to be ^ = {.t|0 < X < \} BXid B = {x\\ < x < 1} 
without loss of generality. 

Fig. 1(a) shows the probability of projection onto the 

balanced k = N/2 bipartite states for two, four, and six 
bosons at arbitrary repulsion strengths. For free bosons 
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FIG. 1. (Color online) (a) The probability of finding an equal 
number of particles in each half-ring A and B in the Lieb- 
Liniger ground state increases monotonically with repulsion 
strength and decreases with particle number, (b) With = 4 
bosons the balanced case fc = 2 is the most probable situation 
for all repulsion strengths. Black tick marks denote the free 
and impenetrable limit values. 



(c = 0) this probability is equal to 
(Fig. 2) and vanishes as ~ y/oJijiN) in the thermody- 
namic limit N ^ oo. This is the expected result because 
each independent particle may be found in either half- 
ring with equal probability. 

As the interparticle repulsion is turned on, however, 
correlations arise between the positions of the bosons and 
the corresponding success probability deviates from the 
free boson value. These quantum correlations give rise 
to fluctuations in the number of particles in each parti- 
tion that result in the probability p{N/2) increasing with 
repulsion strength c. In the limiting case c — > oo of im- 
penetrable bosons (also known as the Tonks-Girardeau 
(TG) limit) this probability reduces to the analogous 
projection probability in a free spinless fermion gas [23]. 
The characteristic function of this probability distribu- 
tion of particle numbers in finite regions of an infinite line 
[24, 25] and a ring [26] are known. We find that a Gaus- 
sian approximation (details are given in the Appendix) 
to the probability distribution p{k) asymptotically yields 
a balanced projection probability of 



p{N/2) 



21n(2A^e'r+i)' 



^> 1. 



(8) 



This probability decays to zero logarithmically for large 
iV — >■ cx) (Fig. 2). It is clear that this probability distri- 



c = 




FIG. 2. (Color online) Balanced projection probabilities for 
the free boson (c = 0) and impenetrable boson (c — >■ oo) cases. 
Our numerical values (dots) are compared to the exact result 
for free bosons (solid line) and the asymptotic result Eq. (8) 
for impenetrable bosons (dashed line). 



bution is symmetric about A^/2 because p{k) = p{N — fc), 
and we further expect it to be unimodal. 

For general k the projection probabilities in the free 
boson limit are equal to 



p{k) 



1 

2N- 



0, 



(9) 



while in the TG limit we expect this distribution to be 
peaked and centered about k = N /2 for a fixed large N . 
Therefore, at these two extreme limits the most probable 
result of the projective measurement is the balanced bi- 
partite state xab{N/2). Since there are no critical values 
of c G (0, oo), we do not expect this trend to change for 
arbitrary finite repulsion strengths and we argue that the 
balanced projection is the most probable outcome for all 
non-negative values of c. We confirm this conjecture for 
= 4 bosons and show in Fig. 1(b) all possible projec- 
tion probabilities for arbitrary repulsion strengths. 



IV. VON NEUMANN ENTROPY AND 
EXTRACTABLE ENTANGLEMENT 

We now calculate the entanglement between the parti- 
tions A and B by calculating the von Neumann entropy 
5a (fc) of region A. The joint density matrix of the pro- 
jected pure state is pab = \xab){xab\ and the relevant 
reduced density matrix is Pa = ^t^b Pab- In the coordi- 
nate basis where |x) = i{}^xn) ■ ■ ■ ^\xi)\0) with ip'^{x) a 
bosonic field creation operator and |0) the vacuum state 
ket, the joint density matrix may be written as 



PAB{k) 



N 



1 



fc)p(fc)|AA|2 



[V] + [Q] 



{V}{Q} 



X I I I I i^-p(x)FQ(x')e*(^^"'-^e-^'' 
'a J a J b Jb 

X |x)(x'| dx^dx^dx^dx^. (10) 
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FIG. 3. (Color online) The von Neumann entropy between 
half-rings in the projected balanced pure state increases 
monotonically with the strength of repulsion. For the case 
of four bosons the balanced projection fc = N jl is more en- 
tangled than the unbalanced cases for all c (inset). Black tick 
marks denote the corresponding values in the impenetrable 
boson limit. 



The signum functions in the ampUtudes i^-p(x) evaluate 
trivially for each case s'gaixBj—XAk) = 1 so that it factor- 
izes according to F-p(x) = F-p^(x^)i^-pB(xB). Thus, the 
integrals over region B may be evaluated independently 
and tracing away the degrees of freedom in B gives the 
reduced density matrix 



E 

{7'}{C} 



A J A 



X e 



|xa)(x^| dxAdx^, (11) 



where the function G-pi^jC) is 



X / FTO(xB)e'^™-"«dxs. (12) 
Jb 

and we have introduced the permuted momentum vectors 
■^•PA = (A-pi, . . . , A-pfc)^ and \vb = (A-pfc+i, . . . , \-pnY ■ 
The eigenvalues of the reduced density matrix may 
be obtained by diagonalizing the associated homogeneous 
Fredholm integral equation 

ii:(xA,x^)0j(x^) dx^ = aj(/)i(x^), (13) 

which has a degenerate kernel 

K{^A,ylA)^ E Gp(x^)e'^-— GQ(x:4)e-*^e-''-. 

(14) 

Since free bosons condense in the ground state, the 
kernel reduces to a constant when there are no interac- 
tions and the only eigenvalue of the resulting reduced 




FIG. 4. (Color online) In the impenetrable boson limit the 
weighted entanglement f (fc) is a maximum for the balanced 
case k = N 12. The extractable entanglement £pp corresponds 
to £{N/2). Connecting lines serve to guide the eye. 



density matrix is unity. Thus, there is no entangle- 
ment when c — 0. For any non-zero contact repul- 
sion between N bosons, however, there are (^) distinct 
ways of choosing the components of the vectors X-pA and 
\qa- From the observation that the plane waves e*^^'^' 
are linearly independent and square-integrable on the 
support A we conclude that the reduced density ma- 
trix Pa has at most (^) non-zero eigenvalues. Hence, 
the reduced density matrix pA{k) is a rank ma- 
trix and we may posit an upper bound for the entan- 
glement entropy extractable from the projected state as 
SA{k) < log2{iV!/[/fc!(iV - fc)!]} EE S'ub(fc) if we assume a 
flat entanglement spectrum. 

The eigenvalue problem (13) can be solved analyti- 
cally by linear algebra methods for small values of N, 
but becomes increasingly cumbersome for large N. This 
difficulty arises because the projected wavefunctions no 
longer have the convenient structure of the ground state 
Bethe ansatz, that is, the projected states are excited 
and the dimensions of the reduced density matrix scale 
exponentially with particle number. A numerically ex- 
act evaluation of the von Neumann entropy is shown in 
Fig. 3 for two, four, and six bosons. For the balanced 
cases we find that more entanglement is present in the 
projected states as the repulsion strength is increased, 
with significant entanglement produced above the scale 
set by c = 1. In the inset of Fig. 3 we show that we 
can extract more entanglement from the balanced pro- 
jection than from any of the unbalanced cases for = 4 
and any c. This fact is supported, but not rigorously 
proven for arbitrary N, by our calculation in the pre- 
vious paragraph of the upper bound for the von Neu- 
mann entropy: Suh{k ^ -A/2) < 5'ub(A^/2) or, equiva- 
lently, rankp^(fc ^ A/2) < rankp^(A^/2). In fact, for all 
the cases we have considered above the projected states 
Xab are not maximally entangled, that is, Sa < '5'ub- 
Nevertheless, because of the preceding observations and 
the symmetry Sa^^ = Sa{N — fc), we conjecture that 
SA{k ^ N/2) < Sa{N/2) for all c. 

So far we have argued that the projection probabilities 
and von Neumann entropies are largest for the balanced 
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FIG. 5. (Color online) The projectively extractable pure 
state entanglement monotonically increases with repulsion 
strength, smoothly transitioning from the free and impene- 
trable boson limits (black ticks). 

projections fc — N/2 at any given value of the repulsion 
strength. Thus, according to the definition (1), the pro- 
jectively extractable pure state entanglement from the 
Lieb-Liniger ground state is given by the weighted en- 
tanglement of the balanced fixed number projections, 

Epp^p{N/2)SA{N/2)=Ek^N/2, Vc>0. (15) 

This statement is trivial for free bosons as we have al- 
ready shown that the von Neumann entropy vanishes in 
all possible projections for any even TV. In the oppo- 
site limit of impenetrable bosons this assertion is verified 
in Fig. 4 where we give numerically exact results for all 
possible cases of up to six bosons. Also apparent in this 
graph is a slower than linear increase in the extractable 
entropy £pp with respect to the boson number TV. 

We also observe that although the probability of suc- 
cessful projection becomes smaller with increasing boson 
number, the von Neumann entropy Sa increases faster 
with N so that the extractable entanglement £pp in- 
creases with both repulsion strength and the number of 
bosons in the ring (Fig. 5). Furthermore, for the few bo- 
son cases we have analyzed here, this increase is mono- 
tonic with respect to both N and c. Hence, more entan- 
glement can be extracted from these projections of the 
Lieb-Liniger gas in the TG limit of impenetrable bosons 
and we may regard the quantity £pp as a probe of both 
quantum correlations and interparticle interactions in the 
ensemble of fixed number projections of the type we con- 
sidered here. 



V. CONCLUDING REMARKS 

We have used the projectively extractable pure state 
entanglement £pp to quantify the entanglement in 
coarse-grained fixed number projections of the Lieb- 
Liniger ground state. This entanglement measure quanti- 
fies the entanglement present in the set of all projection 
outcomes by giving the largest von Neumann entropy 



weighted by the projection probability of a particular 
measurement result. In our case, this maximum corre- 
sponds to the most probable and most entangled projec- 
tion having an equal number of particles in each half of 
the ring. We have seen in our numerically exact few parti- 
cle results indications that the extractable entanglement 
increases monotonically with the strength of repulsion 
and with particle number. We have also observed that 
significant amounts of entanglement can be extracted by 
this projection procedure only in the strongly repulsive 
regime c ^ 1. This increase and subsequent saturation 
of entanglement with interaction strength c has also been 
previously observed numerically in a few-particle Lieb- 
Liniger gas under the framework of single-particle par- 
tioning of the ground state [11]. 

As the impenetrable TG limit corresponds to free 
fermions, our work reveals a fundamental difference in 
the origin of the entanglement present in bosonic and 
fermionic systems. When fluctuations in partition occu- 
pancies are nullified by appropriate measurements, the 
entanglement of free-bosonic systems vanish, while for 
free-fermionic systems they remain non-zero. In the in- 
terpolation between the two cases, as is evident from 
Figs. 3 and 5, the free-bosonic behavior seems to last for 
c nearly up to unity, before a rapid rise in entanglement 
is seen. A qualitative explanation of the difference found 
numerically in this paper is the fact that all bosons are 
exactly in the same state in the c = limit, so that no un- 
certainty remains as to their state in each partition once 
their number in the partition is determined. On the other 
hand, for free fermions, all states up to the Fermi level 
will be occupied resulting in uncertainties (entropies) of 
the state in each partition even when the exact number 
of fermions in the partition is known. We expect these 
features to be preserved at larger particle numbers TV and 
ultimately present itself as an essential difference between 
a Bose-Einstein condensate and a Fermi gas from the 
perspective of quantum information — a suggestion that 
is supported by our results and the observation that the 
Lieb-Liniger model has no density-controlled phase tran- 
sition at any finite positive density. 

It is important to remark that the projected states we 
have described here are not eigenstates of the original 
Lieb-Liniger hamiltonian (2) and will therefore evolve 
non-trivially in time. The results we present here are 
therefore only valid immediately after measurement while 
the spatial partitioning of particles is meaningful; At ar- 
bitrary times after measurement the number of particles 
in each partition will display quantum fluctuations. The 
time evolution of these projected states and the quan- 
tification of the entanglement in them are left as open 
problems for future work. 
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The characteristic function may therefore be approxi- 
mated by the N x N Toeplitz determinant when N ^ 1: 



/(Q)«det(ljv-(l-e*")r). 



(20) 



The matrix T is identical to the single-particle correlation 
matrix of free fermions on an infinite one-dimensional 
lattice upon making the replacement ni/L — )• kpa with 
kp the Fermi wavevector and a the lattice spacing [25] . 

For \a\ < n the Fisher-Hartwig formulas [27, 28] may 
be used to obtain the asymptotic result [25, 26] 



Appendix 

We restore length units in this Appendix. 

We seek an asymptotic approximation to the charac- 
teristic function f{a) = (^e"^ lo i'H'^M^) dx^^ ^^le prob- 
ability distribution p{k) of finding k impenetrable bosons 
in an arc ^ of a ring of circumference L. In the thermo- 
dynamic limit N ^ 00 with finite particle density N/L, 
the characteristic function is equal to the Fredholm de- 
terminant [21] 



/(a)=det(i-(l 



(16) 



where the linear integral operator V acts on the interval 
[—q, q] with q = {N — 1)'jt/L and possesses the kernel 

Let us make a discrete approximation to the integral 

V[F{fi)]{X) by transforming A A„j = (2m — l)7r/L 
and iJ ^ = (2n. — l)7r/i to obtain 



/ V{K,fj,)Fi^)dn^J2^ijF 

J-q 



27r 



N + 1 



Here, the elements Tij of the matrix T are 



sm TTU 



TT{i - j) 



(18) 
(19) 



..../.^^(,^,/(,,))^(, (21) 

^ ' [2Nsm{n£/L)]"^^''' ^ ^ ' 

where G{z) is the Barnes G-function defined functionally 
through G{l+z) = T{z)G{z) and G(l) = 1, with T{z) the 
usual gamma function. Using a small a -C tt expansion 
[25] gives 



, , iNi ln[2iVe^+i sin(7r^/L)l , 
In /(a) ^ — a - ^ ^ 



-l-Oa'^ 



< a <C TT, (22) 



where 7 is Euler's constant. We extract the first two 
cumulants of the probability distribution p{k) from this 
cxprc;ssion and make a Gaussian approximation for the 
case of interest £ = L/2 about the central peak k = N/2 
to obtain 



p{k) 



(23) 



with variance cr^ 



' ln[2iV, 



,7+11 



Physically, the vari- 



ance 0"^ is the fluctuation of particle number in the half- 
ring £ = L/2 about the mean value N/2. The probabil- 
ity of finding exactly A^/2 impenetrable bosons (or free 
fermions) in a half-ring is therefore 



p{N/2) 



21n(2iVe7+i)' 



iV» 1. 



(24) 
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